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We derive exact analytical solutions of the GOY shell model of turbulence. In the 
. absence of forcing and viscosity we obtain closed form solutions in terms of Jacobi el- 

o , 

\ liptic functions. With three shells the model is integrable. In the case of many shells, 

' we derive exact recursion relations for the amplitudes of the Jacobi functions relat- 



=3 



ing the different shells and we obtain a Kolmogorov solution in the limit of infinitely 
many shells. For the special case of six and nine shells, these recursions relations are 
Q ' solved giving specific analytic solutions. Some of these solutions are stable whereas 

^ ■ others are unstable. All our predictions are substantiated by numerical simulations 

of the GOY shell model. From these simulations we also identify cases where the 
models exhibits transitions to chaotic states lying on strange attractors or ergodic 

(N 

J> ' energy surfaces. 

m ■ 
(N 

^ ; I. INTRODUCTION 

in 
o . 

^ ■ The behavior of turbulent systems is fairly well understood from numerical simulations 

of the Navier-Stokes equations but a deep analytical understanding of fully developed tur- 

' n 

^ ' bulence based on these equations is still lacking (see for instance 7] and references therein). 

H : 

^ ■ A fully developed turbulent state is obtained when the Reynolds number of the flow grows 
towards infinity which for instance can be realized in the limit of vanishing viscosity. This 
limit of the Navier-Stokes equations is known to be singular and when the viscosity is identi- 
cally zero the system is usually termed the Euler equations (without external forcing) which 
are believed to exhibit a singularity in a finite time. 

Direct numerical simulations of the Navier-Stokes equations at high Reynolds numbers 



are generally very cumbersome and in the last decades s^ 



lell models of turbulence have 
23|. Indeed, using shell model 



become strong alternatives to direct simulations ,lio| 
approaches one easily simulates Reynolds numbers of the order from 10^° up to 10^^. Shells 
models provide excellent agreement with statistical experimental data, such as structure 



2 



functions, probability distribution functions, etc l|. In contrast shell models generally fail 
to provide knowledge about geometrical structures as these models are too purely resolved 
in Fourier space. 

In this paper we take our starting point in shell models with relatively few shells in the 
Euler case of vanishing viscosity and no external forcing. This might not be a very realistic 
case for a real fluid flow but we are nevertheless able to solve the shell model exactly in this 
limit. In the case of three shells we find the full dynamical solutions in terms of Jacobi elliptic 
functions. We furthermore show that the solutions are stable as is confirmed by numerical 
simulations. In addition, we find that no other solutions are possible meaning that the GOY 
shell model with three shells is here shown to be fully integrable. If we consider this system 
as a low- dimensional version of the Euler equations (possessing the same basic structure and 
symmetries as the real Euler equations), our results lead to the interesting conclusion that 
despite the driving there are no finite time singularities for such a system. 
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23|. 



It is well known that the GOY shell model possesses a period three symmetry 
We therefore attack the model with six and nine shells and find again analytical solutions is 
term of Jacobi elliptic functions. We obtain the interesting results that some of these periodic 
solutions are unstable whereas others are stable. This observation is very much related to 
previous studies on periodic solutions in the GOY model [3], in other turbulent systems 
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20|, |2l| and in the Kuramoto-Shivashinsky equation [4, 8|. Kawahara and Kida 
identified periodic orbits in plane Couette flow UM and van Veen and Kida have looked for 
periodic orbits in a symmetry reduced models of fully developed turbulence pj|,|20|. Toh and 
Itano 19l and Walaffe 2j| have identified periodic structures in channel flows and Eckhardt 

n 

et al has studied pipe flows p]. The Kuramoto-Shivashinsky system is a field equation that 
does not model a turbulent fluid but nevertheless possesses chaotic states developed from 

n 

non-linear mode interactions. Ref. [J] identified a multitude of unstable periodic solution 
in the KS equation. However, Frisch and co-workers found that the phase space of the 
KS equations contains small isolated areas which include stable periodic solutions [8|]. It is 
nevertheless believed that a chaotic attractor of a large dimension in a high-dimentional space 
can be though of as composed of an ordered infinity of unstable periodic orbits ("Hopf's 



last hope", see 



Indeed long unstable periodic orbits have been obtained for the 



GOY model subjected to forcing and viscosity 



13] thus supporting this picture of a strange 



attractor. Our results therefore bridge these observations by finding exact analytical periodic 
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solutions to the unforced GOY model which appear both to be stable and unstable. It is 
important to note that our analytical solutions possess a phase freedom so our results define a 
continuous infinity of exact solutions for the GOY model. Since we shall show (see Appendix 
A) that for infinitely many shells a Kolmogorov solution is obtained, this means that the 
periodic solutions are physically important and constitute a non-trivial part of the phase 
space. 

We also note that one can obtains similar Jacobi type solutions for the Sabra shell model 
jiol although the specific recursion relations for the amplitudes of the Jacobi elliptic functions 
possess different conjugations (see Appendix B). For the Sabra model self-similar quasi- 
soliton solutions have been identified in Ref. 17 1. 

When a small forcing term is added we find that some but not all of the stable solutions 
become unstable. Numerically we observe that the dynamics either becomes chaotic on a 
strange attractor in the 12-dimensional space or quasiperiodic with a regular time evolution. 
The transition can be driven by varying the initial condition in the six'th shell and we show 
numerically examples of such transitions. 



II. SOLUTION OF A THREE-SHELL MODEL 



In this section we shall find solutions to the GOY-model with three shells for the cases 
with no forcing where the phases are constant, and where the phases are variable. This 
model is integrable. Further, we also consider the model with (an imaginary) force, and find 
a solution. The main point in this investigation is that it will serve as a warm up exercise 
for the solution with an arbitrary number of shells given in the next section with further 
details in Appendix A. 



A. Constant phases 

The GOY model is a rough approximation to the Navier-Stokes equations and is formu- 
lated on a discrete set of /c-values, /c„ = r". In term of a complex Fourier mode, m„, of the 
velocity field the model reads 



with suitable boundary conditions [l|, |9|, [lOl, l23|]. / is an external, constant forcing (here 



31] in a reduced version 



on the first mode) and u is the kinematic viscosity. From these equations we have energy 
J2 l^nP as well as hehcity Z](— l)"fcn|MnP (for 3D where 6 = 1 — 1/r) and enstrophy J2 ^nl^nP 
(for 2D where 6 = 1 + conservations 

The mathematics of the three shell case has been studied before 
of 2D MHD with interactions between only three wave vectors. The structure of these 
equations is the same as the three-shell GOY model discussed below. However, the physics 
of the model in ref. [sl is in 2D and is not the same as the GOY model, where we can study 
different dimensions via different conservation laws. 

Initially, we shall for simplicity solve the three-shell model without any force or viscosity 
(i.e. / = z/ = 0) and with constant phases of the velocity functions m„. It is then easy to see 
that the sum of these constant phases must be 7r/2 mod ir. Therefore we replace Ui by iui, 
and take the new ui as well as U2 and to be real. The three-shell equations then become 

^ = ru^it) usit), (2) 
^ = -r6u^{t)us{t), (3) 



and 

-r{^-6)u,{t)u2{t). (4) 



dt 

From these equations we have energy conservation as well as helicity (3D) or enstrophy (2D) 
conservations. However, the three-shell model has two other conservation laws. Combining 
([3]) and (jlj) we easily find, 

u,{t)' = L + ^U2itf. (5) 



Here L is a constant which is determined by the initial conditions. In a similar manner we 
obtain from eqs. ([2]) and ([3]) 

u,{ty = M-^U2{tf, (6) 
where M is again a constant determined by the initial conditions. Thus, 

L = usi0f-^u2{0f, M = u,{0f + ]M0f, (7) 

We note that L + M is the total energy. 

The conservation laws ([5]) and imply that the three shell model with constant phases 
and no forcing is completely integrable (see eq. ([9]) below). Hence there is no chaotic 
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behavior. Later we shall show that even if a constant imaginary forcing is included, the 
model is fully integrable. 
From (131) we have 

-r5 dt. (8) 



By means of eqs. ([5]) and (E]) we can express ui and in terms of U2 and eq. ([8]) then gives 

^"^'^ :-r5t. (9) 



5 "2 



This solution is thus in general an elliptic integral. If we substitute 2J] = 5M = cnx, 
where cn is the Jacobi elliptic function, we obtain 



U2{t) = ^6\u,{0)\^ + \u2m' 



X cn(v^rv/K(0)P + (l-^)K(0)Pt + cn-M )), (10) 

V V\/%i(0)P + k2(0)|V; 

where the modulus of the Jacobi elliptic function is given by 

,2_ (i-^)Af _ K(o)p + i/^K(o)P . . 

L + (1-5)M ^ ^|«3(o)|2 + (i_5)|mi(0)|2" ^ ^ 

Also, cn~^ is the inverse Jacobi function. The Jacobi functions depend on the elliptic 
integrals {k'"^ = 1 — k'^) 

K ^ K{k) = r" , K' ^ K{k'). (12) 

"'o VI - k^ sm"^ 9 

By use of the conservation laws eqs. ([5]) and we can easily find the other u's, 



X sn (v^rv/M0)P + (l-5)K(0)P t + cn~' ( ^— ==— = | 1 , (13) 

V VV'^|«l(0)P+|M2(n^|2 

where we reintroduced the phase in ui, and 



usit) = V'|n3(0)|2 + (l-5)|ni(0)|^ 



X dn ( v^r7|M3(0)|2 + (l-5)|Mi(0)|2 t + cn-^ ( ""'^"^ 1 1 . (14) 

* ^ * /5|mi(0)|2 + |«2(0)P ' ' 



Here we used that the Jacobi elliptic functions satisfy 

sn^ X + cn^ a; = 1, dn^ x + k'^sn^ x = 1, (15) 
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and their derivatives satisfy 

d sn X d cia X d dn x n / \ 

— ; = cn xdn x, — = — sn xdn x, and — ; = —k sn xdn x. (16) 

dx dx dx 

These relations clearly shows the intimate connection between the Jacobi elliptic functions 
and the three-shell model. In the next sections we shall see that this remark can be gener- 
alized to an arbitrary number of shells. 

As to the integrability of three-shell model, one can argue geometrically as follows: Since 
the total energy is conserved, the time evolution is confined to the surface of a sphere in 
R^. Now take any one of the other two conserved quantities eqs. (l5|6l) . This corresponds 
to a cylinder along the Ms-axis with elliptical cross section. The intersection between this 
cylinder and the sphere gives then the line along which the time evolution can proceed. The 
second conserved quantity could limit everything to points, but it does not do so, because 
it is not linearly independent of the other two. So this is then the geometrical origin of the 
solutions in terms of elliptic functions. 



We mention that the elliptic functions can be defined by (see 22l |) 



u 



Mu,k) dt /•! dt 



y/(l-t2)(l-A;2t2)' 7cn(n,fc) ^(1 - (^/2 + p^2) ' 



and u = [ , =, (17) 

Jdn{u,k) ^(1 _t2)(^2 _ ^/2) 

They have product representations and Fourier expansions 

,1/4,-1/2. ^ 1 - 2g2" cos 2x + g^" 27r ^ q"+'/^ sin(2n + l)x 
sniu,k) = 2q^k ^ smx^ ^ _ 2,2„-i 2x + g-- = ^ ? ' 

(18) 

f M o i/4,-i/2,/i/2 ^ l + 2g2»cos2x + g^" 2n ^ g"+V2 ^os(2n + l)x 
cn{u,k) = 2q^k Ik' cosx\{ ^ _ 2,2.-1 2. + g^^-^ = Kfc ^ YT^^^ ' 

(19) 

and 

l + 2g2"-icos2x-fg" vr , 27r ^ cos 2^3; 
'^<^^k) = k n 1 _ 2g2n-i cos 2x + g^-2 = 2K + TT?^' ^^'^ 

Here x = ttu/2K and g = exp(— Trf^'/i^). The constants K and K' can be obtained from 

the hypergeometric function 

K=\ 2F,{\, 1; k') and K' = ^ ,F,{\, 1; 1; 1 - k% (21) 



7 



The hypergeometric function 2-^1(0, 6; c; 2;) is analytic in the plane cut from from z = 1 to 
00. Since k"^ depends on 6 through eq. ffTTl) all results can be analytically continued in the 
cut plane from 2D to 3D. 




V "V 
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FIG. 1: Plot of numerical solutions to a GOY shell model with three shells and 
complex amplitudes in the absence of forcing and viscosity. Plotted are: 

Re ui(t)(red), Re U2(t) (green), Im U3(t)(blue) with the initial conditions: Reui(O) = 0, Imui(O) = 
0, Re U2(0) = 0.5, Im U2(0) = 0, Re U3(0) = 0, Im 113(0) = 0.5/VS. a): 5 = 0.5, b): 6 = 0.1. The 
curves correspond to the exact solutions Eq. ([26]) . 



It should be noted that for small k the Fourier expansions in practice only contain very 
few terms, since q is small. This also applied when k is not so small. For example, if A;^ = |, 
then K = K' and hence q = , so is fast decreasing. However, when k approaches 1, 
a large number of terms are needed in the Fourier transforms because in this case q is not 
small. This follows because q can be expressed in terms of k by the formula 

q = X + 2X' + 15A^ + ISOA^^' + A = i ['^^Z^,]^ - (22) 

In the case L = 0, i.e. ^3(0)^ = (1 — 6)/6 ^2(0)^, the solution ([9]) does not give an 
oscillating solution any longer. Instead one obtains 

«2(t) = ^^2(0) 

^ 2 v/ni(0)^ + n2(0)^Ae-V^^^(^-^) ^^^^ 

V^ni(0)2 + M2(0)V5 + wi(0) + (v/«i(0)2 + U2i0)y5 - m{0)) e-^Vmi-s) h 

In this case k = 1, and obviously an infinite number of Fourier modes are needed to produce 
the exponentially decreasing behavior fl23|) . From eq. fl23l) one can easily obtain ui by use 
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of the conservation law and 



(24) 



y'ixi(o)2 + M2(o)V5 + Mi(o) + i^uiioy + «2(0)V5 - wi(0)) e-V^^d--^) 

Thus, at t ^ oo all the energy has gone to the Ui— mode. 

In three dimensions 1 — 5 is positive, whereas in two dimension this quantity is negative. 
From (dl]) it then follows that fc^ > in 3D, whereas fc^ < in 2D if |m3(0)P > {6-l)\ui{0)\^. 
The Jacobi elliptic functions can be continued to negative k'^, so this does not pose any 
problem. The continued functions are given by 

1 sn(xVl + k^, -/=m) cn{xVTT¥, ) 



dn(x,iA;) = , ^ . (25) 

For k —>■ 0± these functions are continuous. 

The solutions ([T0l)-([T4D simplify if ui{0) = and ^2(0)7'^ = ^3(0)^ From ^ we see 
that in this case k"^ = 1 — 6, and 

Ui{t) = i—j= sn(art, Vl — 5), U2{t) = a cn(art, — 5), U3{t) = dn(art, Vl — 5). 

(26) 

Here a = U2{0). In Fig. [1] we show these exact solutions by direct numerical integration 
of the GOY model with three shells (with specific initial conditions). Note the familiar 
behavior of the Jacobi elliptic functions. Two different values of the parameter 6 are shown 
in order to vary the elliptic functions parameter fc^ = 1 — 5. The solutions are stable and 
the numerical integrations can proceed indefinitely. 

We also mention that if = 0, corresponding to 5 = 1, sn and cn become sin and cos, 
respectively, and dn approaches the constant one. From eq. ( |TT]) we see that fc^ = if 5 = 1. 
This case has been studied (including viscosity) in [isl l. 



B. Solution with variable phases 

We shall now discuss the case of variable phases. In forming the energy or other conserved 
quantities they do not enter directly. However, in this section we shall see that the phases 
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do play a dynamical role. Letting the phases be allowed to vary a new solution emerges, 
which can still be expressed in terms of Jacobian elliptic function, but which has one more 
free parameter than the solutions (|T0l) -( fT4l) . 

The three-shell model with variable phases defined by m„ = iMnle*"^" can be rewritten as 
six equations. 



d\ui\ 



rdt 
and 



\u2\\Us\ Sin(0i + 02 + 03), 



d\u2\ 



d\u3\ 
rdt 



' rdt 

-{1 - 6)\u2\\ui\ Sin(0i + 02 + 03) 



5|Mi||m3| Sin(0i +02 + 0, 



3), 



(27) 



as well as 



rdt 
and I Ma I 



\U2\\U3\ COS(01 + 02 + 03) 



rdt 



-6\Ui\\u3\ COS(01 + 02 + 03) 



rdt 



-(1 - (5)|Mi||m2| COs(0i + 02 + 03). 



(28) 



From the last equations we easily derive 

d{(pl + 02 + 03) 



rdt 



COS(0i + 02 + 





M2I 






Ml 





l^ll 






U2 





\u 


1 


M2 




U3 





(29) 



On the left we can use eqs. ( 1271) in the square bracket to get 



d{4>l + 02 + 03 



dt 



coti 



+ 02 + 



(iln(|ui||u2||^^3| 
Jt 



from which we get 



COS(01 + 02 + 03) 



N 



(30) 



(31) 



P1IF2IIM3I 

where is a constant. We see that if the phases are constants, i.e. = 0, their sum must 
be tt/2 mod vr (i.e. 7r/2, 37r/2,...), as already used in the previous subsection. 

From the cosine in (l3T|l we can compute the sine and insert it in eqs. ( |27|) . By this 
procedure we arrive at 



d\ui\ 

rdt \Ui 
and 



d\u2\ 
rdt 



-6 



1 



P2 



M1PIM2PIM3P — iV^ 



il-6) 



(32) 



rdt ' ' l^sl 

In this equation we have inserted a positive square root, corresponding to a positive sine. 
Since the left hand sides of the eqs. (!32l) can also be negative in some range of t values, it 
follows that for such values of t a negative sign should be chosen for in front of the square 
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roots in eqs. (1321) . corresponding to negative values of the sine. In order not to make the 
notation too clumsy we do not indicate this explicitly. Further from eq. (|27|) we can easily 
derive like eqs. ([5]), ([6]) 

ImsT = ^ + ^^|M2r and = M - ^jual^ (33) 



This allows us to express any |m„| in terms of the other two |m|'s. Using this in eqs. ( I32i) we 
obtain for example for \u2\ the equation of motion, 

d\u2\^ 



2S^J\u2\^{M - \u2\ys){L + |m2|2 (1 - 6)/6) - m 



-d{rt). (34) 



This is again an elliptic type of integral. However, it is slightly more complicated than the 
solutions considered in the previous subsection. 

Inside the square root we have a polynomial of third degree in |M2p- If we introduce the 
zeros of this polynomial, this is the prototype of an integral which can be inverted by the 
use of Theta functions [22;]. However, instead of seeking the most general solution, in the 
following we shall study only a limited range of parameters corresponding to a limited range 
of the initial values of the u's. 

To obtain a relatively simple solution we must remember that the Jacobi elliptic functions 
are not really independent, as is seen from eq. ( JT5l) . Therefore one could attempt to make 
an ansatz in such a manner that inside the square root we have three squares of elliptic 
functions, which requires 

= Ai + Bi sn^(art), \u2\^ = A2 + B2 cn^(art), and \u^\'^ = A3 + ^3 dn^(art). (35) 

With a suitable choice of the A and B constants the square root in eq. fl34l) could then 
produce the product sn(at) cn(at) dn(at) to be matched by the derivative in the numerator, 
dcn^{at) = — 2cn(at) sn(at) dn(at) dt (see eq. ( fT6|) ). 

Instead of proceeding from the integral (IMl) it is actually easier to insert the ansatz (l35l) 
directly in the equations of motion ( l32l) . This gives a number of algebraic equations. When 
we differentiate the functions (l35l) in order to insert them on the left hand sides in (l32l) we 
get a result proportional to -B„ sn(art) cn(art) dn^art) / \un\ with n = 1,2,3 by use of eqs. 
f|T6|) . To match the right hand sides of eqs. (|32|) we therefore need 



niP|u2P|M3P — A^^ oc sn(art)cn(art)dn(art). (36) 
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This means that the following conditions must be satisfied 25| . 



aB, = ^B^B2B^, aB2 = 6 ^B.B^Bs, k'aB^ = {I - 6)^ B^B^B^. (37) 

As before in eq. fl26|) . we take k"^ = 1 — 5. Then the restrictions fl371) are easily solved 

Bi = a^/S, ^2 = a^ Bs = a^/S. (38) 

The remaining conditions then become 

= -A^a^S - A^il - 5)a^/5'' + Ai{l - 5)a^/5, (39) 

which is the requirement that the coefficient of sn^ (art) vanish, and 

= + AiAsa^ + AiA2ay6 + A^a^S, (40) 

which expresses the condition that the coefficient of sn^ {art) vanishes, and 

= A2Asay6-AiA3a^+Asay6-AiA2{l-6)ay6+A2ay6^-Aiay6-Ai{l-6)ay6, (41) 

which expresses the condition that the coefficient of sn^ (art) vanishes. 

The solutions of fl39|) - fHT|) should be such that the square root on the right hand sides of 
eq. (1321) should be real. This is not a trivial requirement. One can study this question in 
general by first solving for A3 in terms of Ai, A2 by use of eq. ( l39l) . and ( 1411) then becomes 
a second order equation which gives A2 in terms of Ai. In this way we obtain 

A, = il-S)iA,-A2/5), (42) 



A2 = 777-^ I a^ + (1 - 6)Ai ± 



2(^1- .;..,^^a4 + (l-5)Ai(^2a2^-3(l-5)AiJ 1, (43) 



as well as 

= Al 6{1 -6)+ Ay {2 -6) + Aia^/5. (44) 

The last restriction can be compared to eq. (13T]) . which expresses the constant in terms 
of the initial values of the phases and the \u\'s. Using eq. fl38|) to find the initial |u|'s and 
comparing this to eq. (IHI) we see that the ansatz (15^ is valid only when cosine of the sum 
of the initial phases is equal to plus/minus one, i.e. 

0i(O) + 02(0) + 03(0) = (mod tt). (45) 
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This does not mean that the phases at later times are restricted in this way. 

The requirements fj42l) . fj43|) and (jHj) do not give acceptable (real) results for all values 
of Ai and a, since we need to ensure that in (jSj) A^^ > and that the square root in (H3l) 
should be real. Let us therefore give a numerical example where the results are good. We 
have found that if 5 = 1 — l/r = 1/2 (3D, with the often used value r = 2) and a = 1 the 
following values are completely consistent (the square root is real!), 

Ai = 0.309, A2 = 0.577, A3 = -0.423, A^^ = 0.7685. (46) 

For initial values of the phases and the |m|'s which are not covered by the ansatz (135|) one 
must go back to the standard elliptic integral flMl) . and treat it in other ways than done 
here. 

The solution ( l35l) can be used to obtain an integral representation for the phases. We 
can insert the cosine from eq. (13T]) in the equations of motion for the phases (l28l) . and we 

get 

ft di 

with similar expressions for the other phases. 

It should also be mentioned that it is possible to take into account the viscosity terms in 
the equations for the phases. The result is 

dt 

with similar expressions for the other phases. Unfortunately we have not been able to solve 
for the modulus \ui\ when viscosity is included. 

In conclusion we have seen that if the phases are allowed to vary, this situation gives 
rize to the new solutions fl35|) . So, although the phases do not appear in the energy and 
other conserved quantities and one might therefore consider them unphysical, this is not 
true because of their dynamical importance. As the phases are allowed to vary continuously 
this indeed outlines an infinity of exact periodic solutions to the GOY models with three 
shells. 
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C. Inclusion of a constant force 

We shall now consider the case where a constant force (/ is taken to be imaginary, / if) 
is included in the first shell, so that ([2]) reads 

^ = ru2{t) u,it) + /, (49) 

whereas the other two equations (I3l) and (jll) are unchanged. For this reason we derive like 
in eq. ([5]) 

usit)' = L + ^U2{tf. (50) 



However, eq. ([6]) is now replaced by 

I - f [ dt' u,{t') = \m-^^ u,{t)\ (51) 

Here the initial conditions again determine the constants M and L. The integral on the left 
hand side can be determined from eq. ([3]) and fl50|) . 



^ 1 du2 du2/dt 
— rdui = — = , . (52) 

Thus, 

Jo V 1 - * «,(o) + ,fW+^i 

Inserting this in eq. ([3]) we obtain by use of eq. fl50|) 

duo /■"2(*) duo 



rS t 



U1U3 J«2(0) 



2 I LS 



«2 + a/^ + y 



(L + 1^uI) \ M-\ul- f / In , 



(54) 

This in principle expresses U2 as a function of t, and ui can then be obtained from (ISTI) . 
whereas M3 can be obtained from fl50|) . However, due to the logarithm in the denominator 
we cannot use elliptic functions to invert eq. (IMll . 

If L = and //r is large and positive, it is possible to obtain an asymptotic expression 
for the integral fl5^ by introducing the variable 

z = M--ul , In^- (55) 
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By a partial integration one then obtains 



-j5il-6) rt = - 



dz 



^{2uy6 + 2f/ir^6il-6)) 



2v^ 



1 + 



(56) 



_2wi/5 + 2//(ry'5(l-5))_ 

Here Zq is obtained from the substitution just by inserting the upper hmit U2{t) instead 
of the integration variable U2 ■ With the assumption U2 <^ / 1 In M2 | we can invert this 
expression and obtain 

U2{t) ^ U2{0) e-v/^) f""'"/^. (57) 

This expression works quite well for large f's, as one can see by solving eq. (H9ll numerically. 
Thus the assumption M2 <^ /| lnM2| is self-consistent. 

For / < we have not been able to find asymptotic results. Numerically one can easily 
see that in this case U2{t) oscillates. In the Fig. [2] we show numerical solutions of the GOY 
model with three shells with a small forcing added. The solutions are surprisingly still stable 
even though energy is now added into the system. We observe the interesting phenomenon 
that the dn-function undergoes a period doubling when the value of the forcing is changed. 
We were however not able to identify a series of period doubling eventually leading into a 
chaotic state. 





FIG. 2: Plot of the Jacobi elliptic functions for three shells with addition of a force, a): / = 
(l + i) ■ 10~2 b): / = (5 + 5i) • with otherwise the same values of the parameters and the 
same initial conditions as in Fig. 1. Note that compared to a) the variation in the dn function has 
undergone a period doubling. In both cases a),b) the solutions are stable. 



15 



III. MANY SHELLS 

In this section we shall extend one of the results obtained in the last section. We start 
by defining the functions 

Joi = sn(art + b), Ja2 = cn(art + b), Ja^ = dn(art + b), 
Ja^ = sn{art + b), Ja^ = cn{art + b), Ja^ = dn(art + b), 
Jaj = sn{art + b), (58) 

where the symbol Ja stands for Jacobi. In general Ja„ = Ja„_3. The Ja's satisfy the 
differential equations 

dJan/dt = e„ Jan+iJcin+2- (59) 

Here 

e„ = +1 if n = 1 (mod 3), e„ = — 1 if n = 2 (mod 3), e„ = —k'^ if n = 3 (mod 3). (60) 

The modulus k'^ is to be determined by the requirement that the functions 

Un = An Jttniart + b) (61) 

should satisfy the shell equations up to some order. Similarly, the constants a, b and A^ 
should be restricted from this requirement. Eq. (|6T1) is of course a generalization of (fTOl) . 
( IT3|) . and ( |T4l) . As an example we show in Fig. [3] two solution to the unforced GOY model 
with six shells. By just changing the initial conditions in the six'th shell, the solution change 
from being quasiperiodic to a chaotic signal. Thus in the Hamiltonian 12- dimensional phase 
space stable periodic and quasiperiodic islands are situated between chaotic areas as is a 
well known scenario for Hamiltonian systems of low dimension. 
The shell model without viscosity and forcing 

c/m* / 5 1 — 5 \ . 

= -ikn \Un+lUn+2 " " Un-lUn+l ^ U„_iU„_2 I . (62) 

is satisfied exactly by the ansatz (l6Ti) . because by direct insertion we get 

LHS = y4*are„ Ja„+i(at) Ja„+2(cti) = (63) 

—ikn \ An+lAn+2JO'n+lJO'n+2 </fln-l ^^n+l 5— ^n-l^ri-2</a,„_i Ja„_2 ) , 
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and using that the Ja's are defined mod 3 in the index we get 

Jttn-i = Jan+2 (in second and third terms), Ja„_2 = Jdn+i (in third term). (64) 

which shows that the Ja's cancel out on both sides of eq. (63). The A's should therefore 
satisfy 

raA^en = -iK \ An+iAn+2 An-iA^+i - ^-^A„„iy4„_2 ] • (65) 



This can be rewritten as a recursion relation 

An+2 



knA. 



A 



(66) 



n+l 
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FIG. 3: Numerical solutions of the unforced GOY model six shells. The initial conditions are: 
a): Re ui(0) = 1.2, Im ui(0) = 0.6, Re U2 = 0.5, Im U2(0) = 0.9, Re U3 = 0.5, Im U3(0) = 
2.0, Re U4(0) = 0.0, Im U4(0) = 0.707, Re U5 = 0.0, Im U5 = .5, Re ue = 0.0, Im ue = 8.0. The 
solution is given by a stable quasiperiodic trajectory, b): Re ui(0) = 1.2, Im ui(0) = 0.6, Re U2 = 
0.5, Im U2(0) = 0.9, Re U3 = 0.5, Im U3(0) = 2.0, Re U4(0) = 0.0, Im U4(0) = 0.707, Re U5 = 
0.0, Im U5 = .5, Re ue = 0.0, Im ue = 0.5. By just changing the initial condition for Im ue, the 
solution changes from a quasiperiodic into a chaotic trajectory. 



We have not succeeded in showing that this relation always has a solution. However, it is 
instructive to compare to the results in the last section. Let us start to solve the recursion 
relation fl66|) . 



A3 



A, 



rA\ 



_ iaAl 
6\A,\^-\A2\ 



(67) 
(68) 
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and 

= rA'MA.r- - + - ^ + 1) - S \A.\') . (69) 

Subsequent coefficients get more complex. The next one is 

^ _ ta{Al + ril-S)Al) ^ tad A* 
^ r^A^ rA2 ' 

Here eqs. ( I68l) and (!69|) should be inserted. 

To gain some insight in the recursion scheme let us start by assuming that there are only 
three shells. Then we should require that A4 = 0, i.e. from 



IA2I = V6 \Ai\. (71) 
Further, from the requirement = one obtains 

= -fcV|Aip + {5^ -5 + 1) lAiH^P - S \A2\\ (72) 



Inserting (17T]1 this gives 



\A 



|2 



(1-5)^, (73) 



which is exactly the same as the previously obtained result ( ITTi) . 

We mention that usually the GOY model is just truncated at some shell number, in this 
case n = 3. However, in dealing with the analytic solutions this is absolutely necessary since 
we need the information contained when the two subsequent shell amplitudes (in this case 
n =4 and 5) are put equal to zero, leading to a value for k"^, which shows that only the 
ffist amplitude is a free parameter. Therefore it is necessary to proceed in the somewhat 
unconventional way followed above, also when we consider an arbitrary but finite number 
of shells, as we shall see in the following. 



With six shells 



26| we need 



A = = + -A4 + -V^, (74) 



and 

iak^Al (1 - 5) 



AjAs = = ^ + ^—^ A,A^. (75) 



Eqs. (1741) and (ffS!) can in principle be used to determine k"^. The algebra is, however, quite 
involved. 
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FIG. 4: The solutions to the GOY model with six shells without forc- 
ing and viscosity for 6 = 0.5 with the following scaling coefficients: 04 = 
0.06527718351453, 05 = 0.6704413705072, = 0.44861490325181 and pa- 

rameter values a = 1,6 = 0,Ai = 1. This leads to the following initial values: 
Re ui(0) = Im ui(0) = Im U2(0) = Im U3(0) = Re U4(0) = Im U4(0) = Im U5(0) = Im U6(0) = 
0.0, Re U2(0) = 0.607820697, Re U3(0) = 1.645222028, Re U5(0) = 0.407507916, Re U6(0) = 
0.73807128. The following shell model variables are plotted: a): Im ui(t), Im U4(t); b): 
Re U2(t), Re U5(t); c): Re U3(t), Re U6(t). 



It is possible to simplify the rather tedious algebra by introducing the following ansatz, 



^4 = Ai, = 05 A2, and Aq = ag A3. 



(76) 



Here we have assigned Ai and A4 the phase tt/2, whereas the other A's are real. In the 
following, we take all the A[s to be real. From eq. fIBSl) we see that 



From eq. 



we similarly get 



Al = {6 - ra^) A{. 



A\ 



(77) 



(78) 
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Next we consider the recursion relation (170|) for Aq which gives 

ai + r{l-6) 6 
"6 = \- -■ (79) 

r-^a^ r 

The conditions for the shell model to stop at six shells are represented by eqs. f TMI) and ( 175ll . 
which lead to 

- = — [d - ra^) — a^, (80) 

and 

(1 — (5) ^4 as = ^ {^r^a\a^ + 1 — 5 + r5 0:4^ , (81) 

respectively. In the last equation we used ]^ from eq. (j78l) . 

From these three equations we can obtain the three a's in terms of the parameters r, b. 
This elimination procedure leads to a sixth order equation for 04, and from the solutions of 
this equation the other a's can then be obtained. 

This procedure leads to extremely messy expressions which are not useful, since we do 
not have analytic solutions of sixth order equations. Instead, it is much simpler to assign 
numerical values for the parameters. Taking r = 2 and 5 = 1/2 (3D) we obtain from eqs. 
( 1791) . ( IHOl) and ( IHTi) by use of NSolve in Mathematica (applied to these three simultaneous 
equations) the following real solutions: 

04 = 0.06527718351453, = 0.6704413705072, ag = 0.44861490325181 

= -0.6788645691146, as = -0.05642693495529, ae = -0.46139658555751, 
a4 = -17.345971520666, as = 1.3995954432919, ae = -1.2098836041347, 
Not acceptable (see remark below eq. (I83l) ) : 

a4 = 0.625261, as = -0.309106, ag = -0.407242 (82) 

From eq. (1771) there are relations between Ai and A2, and using the numbers in (l82l) we get 

Aj/Al = 6-rai = 0.369446, 1.85773, 35.1919, (83) 

respectively. Here the fourth line in eq. (182|) has been ignored, since it does not lead to 
a real ratio A2/A1, as assumed in the above derivation. Thus there only three acceptable 
solutions. 

For fc^ we obtain 

7 2 2 

" 80.2634, -0.33209, 0.390219, (84) 
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FIG. 5: The solutions to the GOY model with six shells without forcing and viscosity for 
5 = 0.5. We apply the following scaling coefficients: 04 = —0.6788645691146, as = 
-0.05642693495529, = -0.46139658555751. Setting a = 10,6 = 0,Ai = 1 give the fol- 

lowing initial conditions: Re ui(0) = Im ui(0) = Im U2(0) = Im U3(0) = Re U4(0) = Im U4(0) = 
Im Us = Im U6 = 0.0, Re U2 = 1.362985377, Re U3 = 7.33683587, Re us = -0.076909087, Re ue = 
-3.385191019 

respectively. 

In Figs. HII5] we show two of the exact solutions to the model GOY with six shells as 
given in (!82l) by the first two set of a„-values. In the first case, shown in Fig. HJ we obtain 
the familiar Jacobi elliptic function with the exact relationships between the amplitude as 
predicted by the theory. Notice that the amplitude of Im ui come to exactly 1 even though 
the initial condition of the mode is set to zero. These solutions are stable. In Fig. O we 
show the the second analytic solution in fl82p . Note we here have set the frequency to a = 10 
in order to obtain a value of the elliptic parameter k less than one. In this case, again, 
the numerical integrations starts out tracing out the exact analytical solutions in terms of 
Jacobi elliptic functions. However, as is clear from Fig. [5] these solutions are unstable since 
after some time they turn into chaotic trajectories. We have experimented numerically with 
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the accuracy of the a„ coefficients and the accuracy of the integration. As is expected, for 
unstable solutions in a chaotic "sea", the exact solution can be sustained numerically longer 
the higher the accuracy. 

The GOY model has been modified by changing the order of the complex conjugates of 
the velocity amplitudes resulting in the Sabra model We have checked that the Jacobi 
elliptic functions are also exact solutions to this model. One can derive recursion relations 
like eqs. fl79l). fl80l) . fl8T]) which have a similar structure but the conjugations and a sign appear 
differently. This is discussed in some details in Appendix B. 

In order to study the stability of the solutions further, we have initiated the integration 
with the exact solutions shown in Fig. Hlbut with the addition of a small force. As is seen 
form the simulation (shown in Fig. [6]), the solution starts out to be stable, looking very 
similar to the exact solution, however in the long run also become unstable and chaotic. 
Nevertheless, we have observed numerically that the existence of this transition depends 
very much on the value of the added forcing /. For some values of /, the solutions remain 
stable although slightly changed away from the exact solutions. For other values of / we 
observe the scenario depicted in Fig. [6l We shall not elaborate further on these transitions 
in the present paper. 




FIG. 6: Solutions to the GOY model with six shells and with the same parameter values as in 
Fig. 4 but with addition of a small forcing Re/ = Im/ = 0.001. Note that the solutions appear to 
be stable for a long time but turn into chaotic trajectories eventually. Plotted are: a) Re ui(t); b): 
Re U4(t). 

It is well known for the GOY shell model that at the parameter value 6 = 1.25 the 
enstrophy H = Y.n=i 

conserved quantity. This point therefore corresponds to 
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two-dimensional turbulence with a forward cascade of enstrophy and an inverse cascade of 
energy We have therefore also looked for exact solutions to the GOY model with six 
shells at this point and by using r = 2 and 6 = 1.25 we obtain as before from eqs. (|79l) . (IHOll 
and (IHTl) by use of NSolve in Mathematica the following real solutions: 

a4 = 0.289785, ag = 2.21803, = 0.613153 

^4 = -0.108771, as = 0.112449, = -0.0517172, 

a4 = 0.0856125, ag = 0.160463, as = 0.302194, 

04 = 0.0982932, = 0.359568, = 0.485351 (85) 

with the following relations between the amplitudes A2 and Ai 

Al/Aj = 0.67043, 1.46754, 1.07878, 1.05341 (86) 

and the elliptic parameter k determined by: 

k^a^/Al = -4.9655, -0.473003, -0.0909193, -0.145639 (87) 

Of these four solutions to the 2D analogy for turbulence we have found numerically that 

number one, three and four on the list are unstable whereas number two is stable. As 

examples of this we show in Fig. [7] some of the solutions, the ones corresponding to the 

second and the fourth in the list. As compared to the 3D case, we observe that the solutions 

in Fig. [Tb become unstable even faster and that they clearly behave very chaotically. 

In Appendix A, we present solutions for the case of nine shells. Furthermore, we derive a 

general recursion relation which allows to solve the coefficients for any number of shells. 

If the number of shell goes to infinity one can see from these recursion relations that the 

_ i_ 

only scaling law solution behaves like an ^ kn^ as expected from the Kolmogorov theory 

Q. 

We end this section by pointing out that the solutions based on the ansatz (176!) are 
special. Presumably the general recursion relations fl66|) for the amplitudes of the Jacobi 
functions can have more general solutions, where the constant phases could also enter the 
recursions. In any case, the solutions that we found represent an infinite set of solutions, 
since the continuous parameters Ai, A2, a, b can vary freely. Alternatively, an analysis of 
eqs. ( I6T!) . (l66l) shows that the free parameters are Mi(0), ^2(0), ^3(0) and Ai. Since these 
quantities in general are complex, this represents an cxD^-continuous set of parameters. 
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FIG. 7: Plot of the Jacobi elliptic functions for six shells for the case of 5 = 1.25 for enstrophy 
conservation corresponding to 2D turbulence for two different sets of a„ values from ()85p (the 
initial conditions are adjusted accordingly), a): a4 = —0.108771, = 0.112449, oq = 

—0.0517172. Here the solutions in terms of elliptic functions are stable, b): 04 = 0.0982932, 05 = 
0.359568, oq = 0.485351. The solutions 'follow' the elliptic solutions for a few periods and the 
clearly become unstable. 

IV. CONCLUSIONS 



We have derived a series of exact analytical solutions to the GOY shell model in the 
absence of forcing and viscosity. In this model a shell couples to nearest and next-nearest 
neighbor shells and the very construction of these couplings directly suggest that the solu- 
tions should be formulated in terms of Jacobi elliptic functions. We have presented closed 
form Jacobi solutions in the case of three shells at which point the GOY model is com- 
pletely integrable and possesses a continuous infinity of stable solutions. We have further 
derived exact recursion relations for the amplitudes of the elliptic functions in case of 6, 
9, 12 etc shells which is related to the period three symmetry of the GOY model. Some 
of these solutions are stable and others are unstable. In the last cases the trajectories lie 
on chaotic energy surfaces. This picture of an infinity of solutions, some either periodic or 
quasi-periodic and others chaotic is common for Hamiltonian systems where the trajectories 
will move on surfaces of constant energy. Our results therefore complement the picture pre- 
sented for the GOY model with forcing and viscosity obtained in ISj where long periodic 
unstable solutions are found embedded in the large dimensional phase space of the dynamics 
of a GOY model with many shells. In this limit of many shells we have found solutions that 
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behave according to the Kolmogorov scahng theory for energy spectra, which shows the 
physical relevance of our periodic solutions. Indeed, we have also observed that adding a 
small force, thus destroying the Hamiltonian property of the system, changes the solutions. 
Some of them are destabilized but some might remain periodic even when a force is added. 
The case of a forced GOY model is closer to a realistic system for turbulence and it is for 
future research to investigate how our solutions might be relevant for a turbulent state when 
a true energy cascade over many scales is present. 
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Appendix A: Solutions with nine shells 

In this appendix we shall discuss the case of nine shells. Following the same procedure 
as in the main text we obtain (the equations for n < 7 are the same as before, but for 
completeness we include the n = 6 eq. (1791) below): 

a4 + r(l-(5) S 
ae = h -. (88) 

r-^a^ r 

a7 = — T- - ™4) + -"4 + ^ r^^' 89 



and 



and 



and 



as well as 



as = — (— r^a4a5 + 1 — 5 + rbdi^— h -as H — - — asa4, (90) 

a-j 5 1 — 5 , . 

a<^ = — V-afi + — — asae, (91) 

r^ag r r^ag 

- ayag = — [0 - ra^) — ayae, (92) 



(1 — 5)aja^ = — ^ (— r^a4as + 1 — 6 + r6a4). (93) 
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These coupled equations can be solved using the simple Mathematica program (5 = .5, r = 2) 

NSolve[{z == (x + l)/{2^y) + 1/4, r == -y{l/2 - 2x)/{2'^z) + l/Ax + 1/2 V^, 

s == ~{-Axy + 1/2 + x)z/{2^r) + 1/Ay + l/2^xy/r, 

t == r/{2^s) + l/Az + l/2^yz/s, l/Art == s/2^(l/2 - 2x) - l/2hr, 

Ijlrs ==t/26(-4xy + l/2 + x)},{x,|/,^,r,s,t}], (94) 

where the notation is 

q;4 = x^ ot^ = y-i cte = = r, = s, 09 = t, (95) 

One obtains 32 solutions, some of which are complex, and some which leads to A\ < 0. Only 
seven of these solutions are acceptable. The reader who is interested in the actual numbers 
should just run the above program. Increased accuracy can be obtained by putting , d at 
the end of the square bracket in (1941) . where d is the number of wanted decimals. 

For the a's we can also obtain in general recursion relations by use of eq. ( l66l) . Defining 
An = OinAi if = 4 mod(3),yl„ = anA2 if = 5 mod(3) and An = anA^ if = 6 mod(3), 
we obtain for n > 4 

fnttn 5 , ^ a;„_ia„_2 

an+2 = \- - "n-i H o > (96) 

r" '-an+i r a„+i 

where we defined ai = 02 = 03 = 1. Here (n = 0, 1, 2, 3, ...) 

e4+3n = 1, e5+3n = "(5 " ^4), h+zn = "(1 " 5 + r6a4 - r^a^a^). (97) 

In addition there are two "stop-equations" : Suppose we stop at the number of shells = 3A^, 
then we must satisfy that A^^^i = and A3n^iA3n^2 = 0, which lead to the two equations 

5 a3N-20:3N = ^g^Ig (S - ra^) - a3N-30i3N-2, (98) 

and 

(1 - 6) a3N^2<y3N^i = ^"^^3 (1 - 5 + r6aA - r^a^a^). (99) 

By means of eqs. flUB]) . flUS]) and (1^^ it is easy to write down the relevant equations for any 
number of shells. The remaining problem is the numerical solution of these equations. We 
are not able to show the existence of acceptable solutions in general, but have checked that 
they exist for six and nine shells. 
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In case with an infinite number of sliells, N oo, the "stop-equations" (pH!) and (p9|) 
sliould of course not be applied so we are left only with eq. (|96|) . If we make the ansatz 
ttn ~ r^", then the only solution is 

/?n = -3 , (100) 

which is fully consistent with the Kolmogorov theory for turbulence. The first term in Eq. 
fl96|) is totally subdominant with this ansatz. 



Appendix B: Solutions to the Sabra model 

In this appendix we shall derive results similar to those in Section III for the Sabra-model 
[igI . This model is given by 



Un = ir 



1-5 



-Un-lUn-2 



(101) 



We can now easily check that if the ansatz ( !6TI) is inserted in this equation, like in the GOY 
model, the Jacobi functions cancel on both sides and we are left with the following recursion 
relation 



A; 



n+l 



4* 



This is the Sabra- version of the analogous GOY model equation fl66l) . We can then proceed 
to solve for the first few A's, 

As = 77^, (103) 



A*2 ' 



which is similar to eq. (1671) . 



A, 



A, 



rA\A\ 



(5iAr-i^2p) 



(104) 



similar to eq. (jHSD, and 

^5 



A\A^ 



r{A\Y\A^\\6\A, 



\A, 



-a^fc^|Ai|^ + (5^-5 + l)|AihA2r + 5|A2r , (105) 



which is similar to eq. ( l69i) . One can of course continue like in Section III to obtain further 
equations from the general recursion relation 



For the case of three shells we should again require that ^4 = 0, which gives 



\A^\^ = 6\A,\\ 



(106) 
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Requiring further that = 0, we obtain 

a^e = {l-6) \A2\^. (107) 

Eqs. (11061) and (I107P are completely the same as the analogous equations obtained in the 
GOY model in Section III. However, in the case of more shells this is certainly not expected 
to continue. 
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